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Abstract. The spin-down of a millisecond pulsar in a 
compact binary system leads to self-energy losses, that 
cause the binary's orbital period to increase (the effect 
being of first post- Newtonian (IPN) order). If the pul- 
sar's period-derivative is not exceedingly small, this effect 
can become measurable in long-term high-precision timing 
measurements of the binary's orbital motion. We use ro- 
tating compressible spheroids to obtain an approximate, 
explicit expression for the orbital period derivative, due 
to spin-down, valid to IPN order. This expression can be 
used to observationally constrain the pulsar's moment of 
inertia, and, combined with other observations, the high- 
density equation of state of compact stars. We apply our 
expression to representative models of millisecond pulsars 
in binary systems and demonstrate the importance of in- 
cluding higher-order rotational effects as well as the im- 
portance of the choice of equation of state, which both 
have been neglected, so far, in the literature. The com- 
puted increase in orbital period is larger, by as much as 
32%, when the higher-order rotational effects are included, 
while it can change by more than a factor of two between 
different plausible equations of state. 



1. Introduction and Motivation 

A large number of known millisecond pulsars are mem- 
bers of binary systems (see Phinney & Kulkarni 1994 and 
Lorimer 1998 for reviews). Nine new millisecond pulsars in 
binary systems have recently been discovered in the glob- 
ular cluster 47 Tucanae (Camilo et al. 2000), and future 
observations are expected to yield an even larger number 
of such systems. Long-term timing observations of several 
relativistic binary pulsars have allowed the measurement 
of secular changes in the orbital period, caused by emis- 
sion of gravitational radiation, with very high accuracy 
(see Kaspi 1999 for a recent review). 

Various causes for secular changes in the orbital pe- 
riod of a binary pulsar have been examined to date. The 
first and traditional such study, in the case of a two-point- 



mass binary, has been based on the use of the quadrupole 
formula for the loss of orbital energy in the form of gravita- 
tional radiation and of the binary's Newtonian (Keplerian) 
orbital motion (Peters and Mathews, 1963). The above 
lowest-order result has been extended (Spyrou and Pa- 
padopoulos 1985) to the general-relativistic IPN approxi- 
mation, with the aid of a similarly generalized quadrupole 
formula (Epstein and Wagoner, 1975) and of the binary's 
IPN orbital motion (Spyrou 1981, 1983; see also Blanchet 
and Schafer, 1989). Of central importance to the above 
calculations is the fact that both members of the binary 
are treated as realistic bodies with finite dimensions, in- 
ternal structure and internal motions, as opposed to the 
usual Newtonian point mass of the standard two-body 
Kepler problem. More specifically, the members are de- 
scribed in terms of the so-called inertial (or gravitational) 
mass, namely, their total mass-energy, and the correspond- 
ing position and velocity three-vectors. These parameters 
are valid to IPN approximation and, initially, they were 
proposed by Spyrou (1977a, b;1978, see also Caporali and 
Spyrou, 1981 for an inclusion of finite-size effects in a rel- 
ativistic parametrized description of the realistic binary) . 
As a consequence, the binary's total inertial mass man- 
ifests itself in the IPN description of its orbital motion 
and, therefore, any of the members's characteristics of in- 
ternal structure and motions can in principle affect its 
orbital motion and hence orbital period. Especially in the 
case of the binary pulsar B1913-I-16, it was proved (Spyrou 
and Papadopoulos, 1985) that, based on the then available 
observational data (Weisberg and Taylor 1984), the rela- 
tive IPN correction to the orbital-period shortening, de- 
duced from the above mentioned generalized quadrupole 
formula, was about five orders of magnitude smaller than 
the relative observational uncertainty. 

Another cause for secular changes in the orbital period 
of a binary or binary-participating compact star, are evo- 
lutionary changes in the compact star's total self-energy 
(Spyrou 1985, 1987; see also Spyrou 1988, 1999 for the 
importance, as a cause of orbital-period changes, of the 
accretion-induced changes of binary-participating pulsar's 
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period of axial rotation and radius) . Especially in the case 
of the binary pulsar B1913-I-16, it was claimed (Spyrou 
1987) that an accuracy improvement in the next decimal 
place of the then observed values (Weisberg and Taylor, 
1984) of the orbital period's rate of change and of the asso- 
ciated uncertainty, would permit orbital-period shortening 
due to self-energy losses to be distinguished from that due 
to gravitational-radiation losses. 

In this paper, we will focus on the secular spin-down 
of a pulsar (due e.g. to magnetic braking), that can cause 
a decrease in its rotational kinetic energy and thus, a cor- 
responding change in the pulsar's gravitational mass. By 
conservation of the orbital angular momentum, this leads 
to an increase in the orbital period of the binary (Jeans 
1924, 1928, Smarr & Blandford 1976, Damour & Taylor 
1991). If the pulsar's period is only a few milliseconds and 
its spin-down rate is of the order of 10~^^s/s (which is 
not uncommon for millisecond pulsars), then the induced 
secular increase in the orbital period may be measurable 
with current technology. Previous studies of this effect 
have used only lowest-order approximations for the self- 
energy change of the pulsar, to obtain order-of-magnitude 
estimates of the effect (see e.g. Smarr & Blandford 1976, 
Will 1981, Damour & Taylor 1991). For example, the star 
was considered to be spherical, i.e. higher order rotational 
corrections to the shape of the star were neglected, and no 
particular equation of state was chosen. Here, we present 
a detailed description of mass-energy loss from a rapidly 
rotating pulsar, deriving an expression for the orbital pe- 
riod change that can be used to accurately constrain the 
pulsar's moment of inertia, from long-term timing mea- 
surements of the orbital period. Our treatment relies on 
the description of rapidly rotating neutron stars as com- 
pressible ellipsoids (Lai, Rasio & Shapiro, 1993, hereafter 
LRS). We apply our results to typical miUisecond pulsars 
in binary systems and derive constraints on their parame- 
ters, for the effect to be observable. Moreover, our results 
are derived consistently to first post-Newtonian order. The 
relation between mass-energy loss and orbital period in- 
crease, originally derived by Jeans (1924), assumed New- 
tonian gravity and isotropic mass-energy loss. The inclu- 
sion of IPN terms in the description of the binary's orbit, 
changes the orbital period at 2PN order. Thus, Jean's re- 
lation is valid to IPN order. Finally, we show that, includ- 
ing spin-orbit coupling in the description of the binary's 
orbit, changes the above relation by a term that can be 
neglected for typical binary systems. 

The outline of the paper is as follows: In the next Sec- 
tion 2, we derive the self-energy loss of a rapidly rotating 
compact star, due to evolutionary spin-down, while in Sec- 
tion 3, we derive the corresponding orbital period change. 
In Sections 4 and 5 we present lowest-order estimates and 
exact numerical results for the orbital period change and 
we conclude with a discussion of our results in Section 6. 



2. Self-energy loss due to spin-down 

We will derive a relation between the spin-down rate 
of a pulsar and the resulting mass-energy loss, valid for 
rapidly rotating stars. The structure of rotating stars is 
approximated by compressible spheroids (introduced by 
LRS). These are approximate equilibrium configurations, 
the compressible generalization of Maclaurin spheroids 
(Chandrasekhar 1969), constructed using a variational 
method. The equation of state (EOS) has the usual zero- 
temperature polytropic form 



P = Kp 



l+l/r 



(1) 



where K is the polytropic constant and n is the polytropic 
index. A polytropic index of n = 0.5 corresponds to a 
relatively stiff EOS, while a polytropic index of n = 1 
corresponds to a relatively soft EOS. The rigidly-rotating, 
oblate spheroid has major axis ai, minor axis 03 < ai and 
eccentricity 



= v/1 - (as/a.y 



(2) 



If ffi is the baryonic mass of the spheroid and Rq the ra- 
dius of a nonrotating (spherical) star of the same baryonic 
mass, then the gravitational binding energy W is 
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where, 

g(e) = (l-eY/^: 
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In (g|), the ratio of rotational kinetic energy T to \W\ 
depends only on the eccentricity e and is given by 
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The internal energy U of the spheroid is assumed to satisfy 
the virial relation 



n 

and the total self-energy £ — U 
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The moment of inertia is 
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where k„ is a constant that depends on the equation of 
state (see LRS), with typical values k„ — (1; 0.82; 0.65) 
for n = (0; 0.5; 1). Then, the rotational period of the star 
is defined through 



P = 27rv///2r. 



(9) 
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Thus, through (||), (||) and (||), the total self-energy and 
the rotational period can be expressed entirely in terms of 
e, m, Ro and k„. 

As the star spins down due to magnetic braking, it 
follows a quasi-equilibrium sequence of constant baryonic 
mass (we do not consider significant mass-loss from the 
pulsar). Along this sequence, the total self-energy and ro- 
tational period, for a given equation of state, are then 
functions of the eccentricity only. In order to derive useful 
relations, we employ a series expansion of the structure of 
rotating compressible spheroids, around the nonrotating 
configuration, in powers of the eccentricity e, and retain 
terms of order at most six in e, a quite accurate, indeed, 
approximation for rotational periods of a few milliseconds. 
The total self-energy and angular velocity fl ~ 2Tr/P are 
then written as 
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Inverting the series expansion of Q, one obtains the fol- 
lowing expansion for the eccentricity: 

397n2 + 414n + 225--4 
2800(3 -n)2 
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where 
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(13) 



is a dimensionless form of the angular velocity. Taking 
time derivatives (denoted by an overdot), one can express 
the mass-energy loss £ in terms of the spin-down rate P: 
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where P = 2-k j^ is a dimensionless form of the rota- 
tional period and Jq = |K„mi?Q is the moment of iner- 
tia of a nonrotating star of baryonic mass m and radius 



i?o • With this definition, Rq remains constant during spin- 
down (when changes in the baryonic mass are neglected), 
while the equatorial and polar radii can both change due 
to the changing eccentricity. 

3. Orbital Period Change 

The gravitational mass m of a compact star (entering 
the gravitational equations of motion) can be decomposed 
(Spyrou 1977a) as 



£ 



(15) 



where m is the baryonic rest mass and £ is the Newtonian 
self- energy, i.e 



£^W+U+T 



(16) 



We consider the IPN motion of a binary system, consist- 
ing of two compact objects of gravitational masses irip 
(primary) and rric (companion). The orbital period of the 
binary system of total mass M = nic + rrip and orbital 
energy E = ^Eq (where Eg is the orbital energy per unit 
reduced gravitational mass) is (Blanchet & Schafer, 1989) 
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and fi — TOc?Tip/(?7ic + Trip) is the reduced gravitational 
mass of the system. If the primary loses energy isotropi- 
cally and the orbital angular momentum is conserved, then 
it has been shown by Jeans (1924,1928), assuming mass- 
energy equivalence and Keplerian orbital motion, that the 
induced change in the binary's orbital period is 






(18) 



It is easy to show that under our assumptions, the IPN 
term on the r.h.s. of ( [l7| ) does not contribute to IPN order 
indll). 

We notice that, the pulsar's spin can be included in the 
equations describing the orbital motion, following, e.g. the 
formalism of Gergely, Perjes & Vasiith (1998). In the Ap- 
pendix, we show, that including the spin-orbit coupling, 
the total change in the binary's orbital period is 



Pl^_^rhp 
Ph M 
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(19) 



where rj = fnc/rn-p, ip is the angle substended by the or- 
bital angular momentum L and the pulsar's spin S, a 
and e are the semi-major axis and eccentricity of the or- 
bit and S is the magnitude of the spin. In dl9), only the 
dominant, to IPN order, spin-terms are included. In the 
Appendix, we also show that, for typical millisecond bi- 
nary pulsars, the spin-term in the r.h.s. of (|l9|) is of order 
(i?o/a)^^^ smaller than the first term. Since we will not be 
concerned with binaries that are near coalescence, in the 
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remainder of this paper, we will neglect the spin-term in 
( p^ ) and use only equation (|l^). 

Assuming that the spin-down is mainly due to mag- 
netic braking and, thus, neglecting changes in the baryonic 
mass of the primary, we find 



£ 



(20) 



and through (|4|), ( p^ ) and (|T8|), we derive our final ex- 
pression for the orbital change due to spin-down 
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This expression depends only on the spin and spin-down 
rate of the primary, the total mass M of the binary sys- 
tem and the moment of inertia of a nonrotating star of 
same baryonic mass as the primary. To lowest order in 
the eccentricity, our result agrees with slow-rotation ex- 
pressions derived previously in the literature (see Smarr 
& Blandford 1976, WiU 1981, Damour & Taylor 1991)[], 
and extends those results to rapidly rotating stars. As we 
will show in the remainder of this paper, the higher order 
rotational corrections that we introduce in (E^) are es- 
sential, if one wants to arrive at quantitative constraints 
on the structure of compact objects, using observationally 
derived data of the above effect. 

4. Lowest Order Estimates 

Before applying ( [2l] ) to specific binary systems, it is in- 
structive to use the lowest-order term 







(22) 



for obtaining a qualitative conclusion about the possible 
types of binary systems, in which the orbital change due 
to spin-down can be observationally important. For the 
effect to be measurable, it must be larger than achievable 
timing precision. For example, the precision with which 
the orbital decay Poth/Porh (due to quadrupole gravita- 
tional radiation emission) of the binary pulsar B1913-I-16 
has been measured is 2.3 x lO^^^s^i (Taylor 1992, 1993; 
Kaspi 1999). 

Assuming a radius of i? = 12km, a spin period of 1.6ms 
and assuming that rric ^ rup, we see that the present ef- 
fect exceeds the above precision, if P > 3x 10~^°s/s. In the 
next section, we will show that including the higher order 
corrections in rotation, relaxes this requirement. Typical 
values of P for pulsars with periods of a few milliseconds 



^ Note that, the first two references contain a missprint that 
was corrected in the third reference. 



are IQ-^^s/s <P < IQ-^^s/s (Phinney & Kulkarni 1994, 
Lorimer 1998). Thus, we conclude that in binary systems 
in which the primary is a rapidly rotating millisecond pul- 
sar with periods of a few milliseconds only and spin-down 
rate larger than a few times lO^^^s/s, the orbital change 
due to the pulsar's spin-down becomes, in principle, ob- 
servable. In contrast, the secular increase in the orbital pe- 
riod of the relativistic binary pulsar B1913-I-16 (due to the 
pulsar's spin-down) is several orders of magnitude smaller 
than the secular decrease caused by emission of gravita- 
tional radiation, because of the pulsar's relatively large 
spin period of 59ms (see Spyrou 1987 for a related discus- 
sion - an upper limit to self-energy-loss effects has still not 
been reached by current observational accuracy). Hence, 
in the case of B1913-I-16 and similar relativistic binary 
systems, the influence of self-energy losses on the orbital 
period are still beyond reach of experimental verification. 

5. Results for Typical Binary Pulsars 

We now apply equation (^l|) to different possible millisec- 
ond pulsars in binary systems. We use the accuracy in 
measuring changes in the orbital period of PSR B1913H-16, 
as a measure of the observability of the changes induced 
by spin-down. Thus, we define the ratio 



A = 



4V^b 



2.3 X 10 



-19c 



(23) 



We note that, the measured value of P/P for millisecond 
pulsars is normally corrected by a galactic acceleration 
term. In the case of PSR B1913-I-16, the uncertainty in 
estimating the galactic acceleration term in P/P is 1.9 x 
lO^^^s"^ (Damour and Taylor 1991), which is slithgly less 
than the presently achieved accuracy in measuring P/P 
for this pulsar.. 

In order to illustrate the importance of including the 
higher-order rotational terms, we define 
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Pi 
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n^ 



(24) 



as being the relative difference between the result (21) 
(which includes terms up to O4) and the lowest order esti- 
mate (p2|). We focus on typical millisecond pulsars of mass 
nip — I.4IM0 using two different equations of state that 
span the expected range of possible realistic EOSs: a rela- 
tively stiff EOS, with polytropic index n = 0.5 and a rela- 
tively soft EOS, with polytropic index n = 1.0. For the stiff 
EOS, we choose the nonrotating radius to be Rq — 12km, 
while for the soft EOS, the nonrotating radius is chosen to 
be Rq — 10km. In this way, the corresponding nonrotating 
moment of inertia Iq is, in both cases, in agreement with 
the relativistic moment of inertia computed by Ravenhall 
& Pethick (1994) for various equations of state. Rapidly 
rotating neutron stars may have mass larger than I.4M0, 
if they are spun-up by accretion. However, the neutron 
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Stiff EOS (n = 0.5) 



Soft EOS (n = 1.0) 
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5 X 10"^** 
1 X 10"^* 
5 X 10"^^ 
1 X 10"^^ 
5 X 10"^° 



7.4 X 10"^** 

1.5 X 10~^** 

7.4 X 10"^^ 

1.5 X 10-1^ 
7.4 X 10"^° 



10 
10 
10 
10 
10 



32.4 
6.5 
3.2 
0.7 
0.3 



Table 1. Relative rate of orbital period change (second 
column), corresponding to different values of the pulsar's 
period derivative P. Two different cases for the pulsar's 
period are considered. The quantity A measures the signif- 
icance of the higher-order rotational corrections, while A 
is the ratio of the second column to the accuracy achieved 
in measuring orbital period changes in the binary pulsar 
B1913-f 16. The pulsar's mass is rrip = 1.41Mq, while the 
mass of the companion is rric — Q.2Mq. Results are for a 
stiff equation of state, of polytropic index n = 0.5. 



star's mass enters in the r.h.s. of eq. ( pi| ) only through 
the ratio Iq/M ^ Rq, which is much more sensitive to the 
radius of the star than to its mass. Thus, it is sufficient, 
in the present context, to consider only one representative 
value for the mass. 

The mass of the companion is fixed at rric = 0.2Mq, 
which is a typical value for white-dwarf companions of mil- 
lisecond pulsars in binary systems (our results are weakly 
dependent on the mass of the companion, if the com- 
panion is a low-mass white dwarf). We do not consider 
the case of binary systems with two neutron stars, as, in 
that case, changes in the orbital period are dominated by 
quadrupole gravitational radiation losses (see Spyrou & 
Kokkotas 1994, for an account of self-energy loss effects 
on the gravitational radiation from neutron star binaries) 

Table 1 displays our results for the stiff EOS. For P = 
1.6ms (eccentricity e = 0.517), A = 56.4 for P = 10"^^ 
s/s, while A = 5.6 for P = 10~^^ s/s is also interesting for 
observations. Notice the large relative difference of A = 
32%, which demonstrates the importance of the higher- 
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17.0 
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1.7 
0.3 
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Table 2. Same as Table 1, but for a soft equation of state 
of polytropic index n = 1.0. 



order rotational effects. For P = 3 ms (eccentricity e = 
0.273), A = 32.4 for P = 5 x IQ-^^ s/s and A == 3.2 
for P = 5 X 10-19 s/s, while, for P < 10" ^^ s/s, the 
magnitude of the effect is less than the presently achieved 
observational accuracy. For the larger rotational period, 
the effect of higher-order rotational terms is reduced to 
A = 10%. 

Table 2, displays the corresponding results for the soft 
EOS. For P = 1.6ms (eccentricity e = 0.332), A = 26.1 
for P = IQ-is s/s, while A = 2.6 for P = lO^^^ s/s, 
roughly a factor of two smaller than for the stiff EOS. 
The rotational effects are reduced to A = 19%, For P = 3 
ms (eccentricity e == 0.175), A = 17.0 for P = 5 x lO^^^ 
s/s and A = 1.7 for P = 5 x lO^^^ s/s. The rotational 
effects are now only A = 6%. 

It is evident that the results are sensitive to the chosen 
equation of state. The present uncertainty in the high- 
density equation of state of neutron star matter (which 
could be as stiff as a n = 0.5 polytrope or as soft as a 
n — 1.0 polytrope) introduces, in the computed orbital 
period change, an unknown factor, which could be as large 
as a factor of two. A measurement of the orbital change 
due to the pulsar's spin-down and a knowledge (or guess) 
of it's baryonic mass could yield a value for (or constrain) 
the moment of inertia /q of nonrotating neutron stars of 
same mass m-p. Combined with a measurement (or ed- 
ucated guess) of the pulsar's mass, this information can 
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then constrain the very high-density equation of state of 
compact stars. 

6. Discussion 

We have extended previous, approximate, calculations of 
the orbital period change in binary pulsars (due to the 
pulsar's spin-down), by including higher-order rotational 
effects and the influence of the choice of the high-density 
equation of state. We have shown that the above higher- 
order terms, as well as the choice of the equation of state 
are significant, if quantitative constraints are to be drawn 
for future observations of this effect. We also showed, that 
the inclusion of the pulsar's spin in the description of 
the binary's orbit, does not contribute significantly to the 
orbital period change. Our treatment relies on the de- 
scription of rotating stars as compressible perfect fluid 
spheroids (ignoring the possible influence of a solid crust), 
in which one can express both the change in spin-period 
and the change in radius, during spin-down, through the 
change in the eccentricity of the spheroid. Binary pulsars, 
that are interesting for the above effect to be measurable, 
must be rapidly rotating, with a rotational period of only a 
few ms and, in addition, with a period derivative equal to 
or larger than a few times 10^^°s/s. As Phinney & Kulka- 
rni (1994) point out, a millisecond pulsar like B1937-f21 
(which has P = 1.56ms and P = lO^^^s/s), if it were in 
a binary system, would yield a measurable orbital period 
change. Using our results, i.e. equation (|2l]), one could 
then derive quantitative constraints on the pulsar's mo- 
ment of inertia. More accurate constraints could be de- 
rived, using numerical models of rapidly rotating neutron 
stars that satisfy the general-relativistic field equations, 
for various realistic equations of state. Such a computa- 
tion, however, will make sense only after observational 
measurements of the predicted effect will start becoming 
available. The increasing number of millisecond pulsars 
discovered in binary systems (Camilo et al. 2000) is a pos- 
itive sign that this could happen in the near future. 
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Appendix 

We have shown that ( p8| ) depends only on the Newto- 
nian orbital characteristics (as the IPN orbital terms con- 
tribute at 2PN order). Thus, the pulsar's spin can be in- 
cluded in the description of the binary's orbit, following 
the formalism of Gergely, Perjes & Vasiith (1998), who 
consider the correction to the Newtonian orbit, induced by 
the spin. We define 77 = nic/mp, -0 the angle substended 
by the orbital angular momentum L and the pulsar's spin 
S, a and e are the semi-major axis and eccentricity of the 
pulsar's orbit and L, S the magnitudes of its orbital an- 
gular momentum and spin. Then, the binary's period Pf,, 
orbital energy E and semi-major axis a, are given by 
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and 
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In the case that the spin is zero and the star loses 
mass-energy isotropically (e.g. due to thermal emission), 
it has been shown by Jeans (1924, 1928) that the product 
A = aM remains constant during the orbit's evolution 
(which leads to equation ( p^ in the text). We treat the 
inclusion of spin in the description of the orbital motion as 
a first-order perturbation about the orbital motion with 
zero spin. If so, we can write 
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and £0, Aq are the the values of E and A when spin-orbit 
coupling is ignored. Then, 
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Keeping only first-order terms in the spin (and keeping 
only IPN terms), we find 
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In order to arrive at an order-of-magnitude estimate for 
the spin-term in (p3[), we substitute the lowest-order (in 
spin-period) expressions derived in the main text (using 
S = m), and find 
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Since the spin-term is of order {Ro/a)'^^'^ smaller than the 
isotropic term in (|34|), it can be neglected for typical bi- 
nary pulsar systems considered in the present paper. 



